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Abstract 

Dark Matter direct detection experiments are able to exclude in- 
teresting parameter space regions of particle models which predict an 
important amount of thermal relics. We use recent data to constrain 
the branon model and to compute the region that is favored by CDMS 
measurements. Within this work, we also update present colliders con- 
straints with new studies coming from the LHC. Despite the present 
low luminosity, it is remarkable that for heavy branons, CMS and 
ATLAS measurements are already more constraining than previous 
analyses performed with TEVATRON and LEP data. 
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1 Introduction 

Identifying the nature of Dark Matter (DM) is a central problem in contem- 
porary physics. The nature of DM and how it fits into our current under- 
standing of elementary particles, is not known yet. Numerous astrophysical 
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and cosmological data require the existence of a DM component, that ac- 
counts for about 20% of the energy content of our universe. Although there 
are other possibilities [I], DM is usually assumed to be in the form of sta- 
ble Weakly-Interacting Massive Particles (WIMPs) that naturally freeze-out 
with the right thermal abundance, with a mass of the order of the electro- 
weak (EW) scale. The most studied DM candidate is the neutralino, which 
can be identified as the lightest supersymmetric particle in many supersym- 
metric models. These models own a rich phenomenology due to the large 
number of new particles present in the theory [2]. 

Besides, the existence of extra dimensions are well motivated theoretically 
and DM candidates arise associated to new degrees of freedom. In particular, 
branons have been proposed to explain the missing matter problem within 
the so-called brane- world scenario (BWS) |3]. In this framework, the SM 
fields are forced to live on a three-dimensional hypersurface called brane, 
while gravity propagates on the higher D = 4 + N dimensional bulk space. 
The fundamental scale of gravity is not the Planck scale M p anymore, but a 
new scale Mjj which in this work will be considered arbitrary [3J. 

In this scenario, the existence of extra dimensions generates new fields 
on the brane, giving rise to a Kaluza-Klein (KK) tower of massive gravitons. 
The brane has a finite tension f 4 and its fluctuations are parametrized by the 
so-called branon fields ir a . These branons are the massless Goldstone bosons 
arising from the spontaneous breaking of the exact symmetry existing in 
the case of translational invariance in the bulk space. In the general case 
of translational invariance explicitly broken, the branons are expected to 
be massive fields. When these branons are considered, the coupling of the 
Standard Model (SM) particles to any bulk field is exponentially suppressed 
by a factor exp[—M] (K A 2 /87rf A ], where M^k is the mass of the corresponding 
KK mode of the tower of massive gravitons and A is the cut-off of the effective 
theory that describes the branon phenomenology [HE]. 

If the tension scale / is much smaller than the other new scales so that 
/ 2 C A Mkk-i then the KK modes decouple from the SM particles, so that at 
low energies the only brane-world related particles that must be taken into 
account are branons. In the following we will assume this to be the case and 
accordingly we will deal only with SM particles and branons (3j. 

A could represent the width of brane or any other mechanism that mod- 
ified the short-distance theory to cure the ultraviolet behavior of branons. 
However, for our purposes, A is just a phenomenological parameter. From 
the point of view of the effective theory, A/f parameterizes how strongly 
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(or weakly) coupled quantum brane is, and therefore controls the unknown 
relative importance of tree-level versus loop branon effects. In [3] it was 
shown that the perturbative loop analysis only makes sense for approximately 
A < 47T 1 / 2 / A r ~ 1 / 4 . We will work in such a limit. 

These branons are expected to be nearly massless and weakly interacting 
at low energies. In general, translational invariance in the extra dimensions is 
not necessarily an exact symmetry, so that explicit symmetry breaking leads 
to a branon mass M. Brane fluctuations could be candidates for the cosmo- 
logical DM and they could also make up the galactic halo and explain the 
local dynamics. Therefore, they could be detected by DM search experiments 

Several experiments have been developed to detect DM directly and in- 
directly. The direct DM detection experiments are designed to observe the 
elastic scattering of DM particles with nuclei while indirect DM searches 
may detect the DM annihilation productions such as protons, antiprotons, 
electrons, positrons, neutrinos and gamma rays. Complementary to these, 
collider DM searches are performed at experiments like CERN (European 
Organization for Nuclear Research) LHC (Large Hadron Collider) [HE]. 

The outline of the paper goes as follows: first we shall discuss briefly the 
model in Section [2j Section [3J is dedicated to direct DM search experiments, 
reporting some DM candidate events and their relation with branons. Sec- 
tion H] is devoted to the analysis of branons related to collider results and 
new LHC constraints. Finally section presents our conclusions. 

2 Flexible Brane- Worlds 

The presence of branes in extra dimensional models have been studied from 
many different points of view (read, for instance, [H]). We will consider a 
single-brane model in large extra dimensions, where the four- dimensional 
space-time M4 is embedded in a D-dimensional bulk space. This space is 
assumed to have the form Mp = M4 x B where the B homogeneous space 
is an N-dimensional compact manifold, such that D = 4 + N. The brane 
lies along M4 and its contribution to the bulk gravitational field is neglected 

mum- 

Branons couple to the conserved energy-momentum tensor of the Stan- 
dard Model evaluated in the background metric Tg V M . The lowest order ef- 
fective Lagrangian is given by [3]: 
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C Br = ^d^d^-^MW 

+ w - M2 ^ a ^) n v M . (i) 

Branons are stable and difficult to detect because they always interact 
by pairs and their interactions are suppressed by the tension scale /, which 
implies that they can be weakly interacting. Also they are expected to be 
massive so that their freeze-out temperature could be relatively high and 
then their relic abundances could be cosmologically important. This implies 
that branons are natural candidates to DM in a scenario where / <C [3]. 

The thermal relic branon abundance is calculated in [3] where it has been 
considered the relativistic (hot) and nonrelativistic (cold) cases at decoupling. 
The allowed region for hot branons masses are much smaller than those in 
neutrino dark matter models so that they decouple much earlier than neutri- 
nos; therefore, hot branons are disfavored. Branons could be responsible for 
the observed cosmological DM density provided Qsrh 2 = 0.110 ± 0.006 [llj. 



3 DM Direct detection experiments 

Several DM direct detection experiments have presented recently some re- 
sults. Their reported limits are compatible with the branon scenario. 

Assuming that the DM halo of the Milky Way is composed of branons, 
its flux on the Earth is of order 10 5 (100 GeV/M ) cm s , and could be 
sufficiently large to be measured in direct detection experiments such as 
DAMA, XENON100, CoGeNT, CDMS II or EDELWEISS-II. These experi- 
ments measure the rate R, and energies En of the nuclear recoils |12| . 

The differential counting rate for a nucleus with mass is 

dR po f°° do B 



dE d m N M 



^ vf(v)^(v,E R )dv, (2) 



where p is the local branon density, {dasr-N / dE^jiy, E B ) is the differential 
cross-section for the Branon- nucleus elastic scattering and f(v) is the Branon 
speed distribution in the detector frame normalized to unity. 

The relative speed of a dark matter particle is of order 100 km~ 1 s~ 1 , 
so the elastic scattering is non-relativistic. Then the recoil energy of the 
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nucleoli in terms of the scattering angle in the center of mass frame 6*, and 
the branon-nucleus reduced mass Hn = Mm^jiM + m^), is given by 



Er = P> 2 (1 -cos6*) 
R rn N 

The lower limit of the integration over the WIMP speed is given in terms 
of the minimum branon speed which can cause a recoil of energy Er and is 
given by v m i n = (mNER/2fi 2 N ) 1 ^ 2 . The upper limit is infinite; however the 
local escape speed v esc is the maximum speed in the Galactic rest frame for 
WIMPs which are gravitationally bound to the Milky Way. The standard 
value for this scape velocity is v esc = 650 kms^ 1 [12J. 

Integrating the differential event rate over all the possible recoil energies, 
it is possible to find the total event rate of branon collisions with matter per 
kilogram per day, R. The smallest recoil energy that the detector is capable 
of measuring is called threshold energy, Et- In terms of this threshold energy, 
Et, the total event rate R has the form |12] : 

R = dE R ^- r vf (v) (v, E R ) dv. (4) 

J Et m NM J Vin . n dE R 

The branon-nucleus differential cross section contains the particle physics 
inputs and depends on the branon-quark interaction given by ([T|). For a 
general DM candidate, its nucleus cross section is separated into a spin- 
independent (SI) scalar contribution and a spin-dependent (SD) one |12) : 

Here the form factors Fsi(Er) and Fsd{Er) account for the coherence loss 
which leads to a suppression in the event rate for heavy WIMPs on nucleons 
and includes the dependence on the momentum transfer q = ^IvtinEr. Uq 
and <7q d are the spin-independent and spin-dependent cross sections respec- 
tively, at zero momentum transfer. These quantities still depend on nuclear 
structure through isospin content; that is, the number of protons vs. neutrons 



In the branon case, the entire interaction is SI and can be written, in 
general, as [13] 



[Zf P + (A- Z)f n } 2 
a = - n 72 a p , (6) 



fp I 1 DM 



V 



5 



with A the atomic mass number, Z the charge of the nucleus, f P:H the SI 
DM couplings to proton and neutron respectively. [iDMp is the reduced DM- 
proton mass, and a p 1 the SI cross section for scattering of DM on a proton. In 
particular, branons do not violate isospin symmetry if we neglect the differ- 
ence in mass of protons and neutrons. Therefore, within this approximation: 
f p = f n , and yU = fiBm — l^Brp- Indeed, the branon-nucleon cross section a n 
can be written as [3] 

SI _ _ 9M 2 m n V 

p ~ n ~ 64tt/8 ' {<} 

where m n is the nucleon mass. Recently, direct search experiments have 
reported possible candidate events for DM. The annual modulation signature 
found either by former DAMA/NAI and DAM A/LIBRA detector localized 
in Gran Sasso National Laboratory points out to a light WIMP [14J. A 
similar conclusion can be obtained from the CoGeNT observations. This 
ultralow-noise germanium detector operated deep underground in Soudan 
Underground Laboratory has found some events consistent with a WIMP of 
mass 7 to 11 GeV |15| . However, these measurements are in clear tension 
with exclusion limits obtained by other experiments located in the same 
laboratories such as XENON100 (a liquid xenon detector at the Gran Sasso 
National Laboratory) [16], or CDMS II (a germanium and silicon detector at 
the the Soudan Underground Laboratory) |17| . 

The CDMS Collaboration have reported two candidate events for DM 
that are consistent with heavier WIMPs and with the present constraints 
|17) . Indeed, the analysis was performed on data taken during four periods 
between July 2007 and September 2008, maximizing the expected sensitivity 
for a 60 GeV WIMP. Their spectrum-averaged equivalent exposure for a 
WIMP of this mass is 194.1 kg-days. Two events in the WIMP acceptance 
region were observed at recoil energies of 12.3 keV and 15.5 keV. 

Neutrons with energies of several MeV can generate single-scatter nuclear 
recoils that are indistinguishable from possible DM interactions. Also an 
approximation made during the ionization pulse reconstruction degrades the 
timing-cut rejection for a small fraction of surface events with ionization 
energy below 6 keV. Therefore, the probability to have observed two or more 
surface events in this exposure is estimated to 20%. Including the neutron 
background, the candidate events have a probability 23% of not being due 
to a DM signal [T7] . 

The results of the analysis for branon-nucleon cross section a n in terms of 
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Experiment 


^(TeV) 


^(pb- 1 ) 


/o(GeV) 


Mo(GeV) 


HERA 1 


0.3 


110 


16 


152 


Tevatron-II 1 


2.0 


10 3 


256 


902 


Tevatron-II 2 


2.0 


10 3 


240 


952 


LEP-II 2 


0.2 


600 


180 


103 


LHC 1 


7 


3.4 x 10" 4 


189 


3240 


LHC 1 


7 


11.7 x 10~ 3 


236 


3240 


LHC 2 


7 


1.21 x 10" 3 


152 


3390 


LHC 1 


14 


10 5 


1075 


6481 


LHC 2 


14 


10 5 


797 


6781 



Table 1: Limits from direct branon searches at colliders (results at the 95 % c.L). Upper 
indices ' denote mono jet and single photon channels respectively. Current data [20l 
I2T1 and prospects for the LHC are compared with present constraints from LEP [7] , 
HERA and Tevatron [S]. y/s is the center of mass energy of the total process; C is the 
total integrated luminosity; fo is the bound on the brane tension scale for one massless 
branon (JV = 1) and Mq is the limit on the branon mass for small tension /—>•[). It is 
important to note that the effective approach taking in this analysis that allows to write 
Lagrangian {TJ is not valid for energy scales A > 47T 1 / 2 / N^ 1 / 4 [5j. Therefore, the Mq 
value cannot be trusted. In any case, we are providing this number since it is the most 
simple way to characterize the sensitivity of the analysis for heavy branons. 

the branon mass and the limits of direct DM search experiments are shown in 
Figure [TJ For reference, lines of constant / with 50 GeV separation are shown. 
The area on the left of the Qsrh 2 = 0.122—0.098 curves is excluded by branon 
overproduction, but the right portion is compatible with observations. Such 
region corresponds to / > 120 GeV and M > 40 GeV. The dashed yellow 
region is favoured by the 2 CDMS II events at 12.3 keV and 15.5 keV at the 
90 % confidence level. 

4 Colliders Searches at the LHC 

If the new physics able to explain the DM puzzle is related to the electroweak 
scale, it may be accessible at the LHC and in the new generation of collider 
experiments [TJ [8] • 

The most important process for branon production in a proton-proton 
collider, as the LHC, is given by the gluon fusion giving a gluon and a bra- 
non pair; as well as a quark-gluon interaction giving a quark and a branon 
pair. In these cases, the expected experimental signal is one monojet J and 
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M(GeV) 

Figure 1: Elastic branon-nucleon cross section a n in terms of the branon 
mass. The two thick lines correspond to the flsrh 2 = 0.122 — 0.098 curve 
for cold branons with N = 1 (left) and N = 7 (right). The shaded areas on 
the left are the collider exclusion regions [6j [8], also for N=l, 7. The solid 
lines correspond to the current limits on the spin-independent cross section 
from direct detection experiments XENON 100 [16], EDELWEISS-II [TS] and 
CDMS \17\ . The striped area is favored by CDMS measurements [17] and 
the dotted lines are sensitivity prospects for XENON 100 and super CDMS. 
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missing energy and momentum which can easily be identified. An additional 
interesting process is the quark- ant iquark annihilation giving a photon and 
a branon pair. For this process, the signature is one single photon and miss- 
ing energy and momentum. All the branons are assumed to be degenerated 
with a common mass M and the quarks are considered massless. Top quark 
production has been neglected due to its large mass. 

The cross-section of the subprocess gg —> gun is given by [5] 



d(j(gg -> gnn) a s N(k 2 - AM 2 ) 2 I ANP 



dk 2 dt 409Q0f 8 7r 2 sHu V k 2 

(s 4 + t 4 + m 4 - k s + 6k\s 2 + t 2 + u 2 ) - Ak 2 {s 3 + t 3 + u 3 )), (8) 

where s = (pi + P2) 2 , t = (pi — q) 2 , u = (p 2 — q) 2 and k 2 = {k\ + /c 2 ) 2 - Pi an d 
P2 are the initial gluon four-momenta, q the final gluon four-momentum and 
k = k\ + k2 the total branon four-momentum. Therefore the contribution to 
the total cross section for the pp —> girir reaction coming from this subprocess 
is given by 

v gg (pv -> Qw) = dx dyg(y\ s)g(x; s) 

%ax Al 2 f tmax i + d(J (.39 -»■ gw) , v 

Here g(x; s) is the gluon distribution function of the proton, x and y are the 
fractions of the protons energy carried by the initial gluons. The different 
limits of the integrals can be written in terms of the cuts used to define the 
total cross-section. For example, in order to be able to detect clearly the 
monojet, a minimal value for its transverse energy Et and a pseudorapidity 
range given by r] min and r) ma x must be imposed. Then the limits fc^ in = AM 2 , 
Kulx = s(l - 2E T /y/s) and t min{max) = -(s - k 2 )[l + tanh (r) min ( max ))}/2 are 
obtained. Also we have x m i n = s m i n /s and y m in = Xmin/x where s is the 
total center of mass energy squared of the process and 



2E1 + 4M 2 + 2E T JE 2 +4M 2 . (10) 
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For the qg — > qirir subprocess the cross-section is given by [5] 



da(qg — > qim) 



dk 2 dt 




with p\ and p 2 being the quark and the gluon four-momenta respectively, 
q the final state quark four-momentum and k\ and k 2 the branon four- 
momenta. The Mandelstam variables are defined as in previous cases. 
The total cross section for the the reaction pp — > qirir is then 



In this equation x and y are the fractions of the two protons energy carried 
by the subprocess quark and gluon. The different limits of the integrals can 
be written as in the previous case in terms of the the minimal transverse 
energy of the quark (monojet) Ej>- 

Considering all the above equations, it is possible to compute the total 
cross section a(pp — > Jim) in terms of the cut in the jet transverse energy 



On the other hand, to analyze the single photon channel, we need the 
cross-section of the subprocess qq — > '-firir, that was computed in [8]: 




(12) 



da(qq — > jim) 



dk 2 dt 
Q 2 q aN{k 2 -4M 2 ) 2 



1 - 



4M 2 



(sk 2 + Atu){2sk 2 + t 2 + u 2 ). (13) 



184320/% 2 s 3 tw 



And this is the only leading contribution: 




(14) 
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By using all these cross sections, it is possible to compute the expected 
number of events in these channels produced at the LHC in terms of the brane 
tension scale /, the branon mass M, and the number of branons N. The main 
source of uncertainty is coming from the parton distribution function that 
we have taken from |19] . 

Previous constraints on the branon model parameter for tree-level pro- 
cesses from other colliders are summarized in Tabled] |H1 [SJ. In this Table, 
the present restrictions coming from HERA, Tevatron and LEP-II are com- 
pared with the present LHC bounds running at a center of mass energy 
(c.m.e.) of 7 TeV and the prospects for the LHC running at 14 TeV c.m.e. 
with full luminosity. For the single photon channel, ATLAS has published 
two different analysis with a total integrated luminosity of 0.68 nb~ l and 
0.53 rib' 1 [20]. The results are not well understood within present hadron 
interaction models, but we can deduce the constraints showed in Table [T] by 
assuming a conservative approach. In any case, it is much more interesting 
the analysis presented by CMS about jet and missing transverse energy, that 
is the most constraining one for branon phenomenology due to its high lumi- 
nosity: 11.7 nb~ l |21j . In order to complete the analysis, we have included 
on Table [TJ the analogous study by ATLAS which is not competitive due to 
the low luminosity of 0.34 nb~ l used in the analysis |22j. 

In addition, it has been shown that branon loops introduce new cou- 
plings among standard model particle that can be described by an effective 
Lagrangian. The most relevant terms of this effective Lagrangian are [5]: 

/VA 4 

£ ^ " { 2T - T " + T ^ ■ (15) 

The A parameter appears when dealing with branon radiative corrections 
since the Lagrangian (CQ) is not renormalizable. This parameter is the cut- 
off which limits the validity of the effective description of branon and SM 
dynamics. A one-loop calculation with the new effective four-fermion ver- 
tices coming from (fT5|) is equivalent to a two-loop computation with the 
Lagrangian in (0Q), and it also allows to obtain the contribution of branons 
to the anomalous magnetic moment of the muon |5]: 

5j iVA 6 

5a » " Tl4W"F (16) 

where N is the number of branon species. 
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Experiment 


(TeV) 


C (pb- 1 ) 


/ 2 /(iV 1 / 4 A) (GeV) 


HERA C 


0.3 


117 


52 


Tevatron a ' b 


1.8 


127 


69 


LEP a 


0.2 


700 


59 


LEP 6 


0.2 


700 


75 


LHC 6 


14 


10 5 


383 



Table 2: Limits from virtual branon searches at colliders (results at the 95 % c.L). The 
indices a ' b ' c denote the two-photon, e + e~ and e + p (e~p) channels respectively [5]. Present 
constraints from HERA, LEP and Tevatron are compared with the LHC prospects. The 
first two columns are the same as in Table [IJ and the third one corresponds to the lower 
bound on f 2 /(N^ 4 A). In this case, current LHC constraints are not competitive. 

The most relevant branon loops that could have compatible effects with 
SM phenomenology could be the four-fermion interactions or the fermion pair 
annihilation into two gauge bosons [5]. Bounds on the parameter combination 
f 2 / (AiV 1 / 4 ) can be established by considering current data. Present results 
are shown in Table [21 where it is also possible to find the prospects for LHC 
[5]. However, in contrast to branon direct production, current measurements 
by CMS or ATLAS detectors are not able to improve the limits obtained 
with data coming from HERA j23], Tevatron |25] and LEP [26J. 

The preferred parameter region for branon physics is given by [5] 

6 -° GeV £ nCa* ~ 2 - 2 GeV ( 95%cl ) ( 1T ) 

As a result of this relation and the limits shown in Table [21 the first branon 
signals at colliders would be associated to radiative corrections j5]. 

5 Conclusions 

Brane fluctuations are DM candidates in brane world models with low ten- 
sion. We have studied the current situation of branons in relation with direct 
dark matter detection experiments and recent LHC measurements. As it is 
possible to observe in Figure [21 the results show that, in a certain range of the 
parameters for the brane tension / and branon mass M, the relic abundance 
could explain the missing mass problem and that such parameter regions 
would be compatible with DM candidate events currently observed by di- 
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2500 
2000 
1500 

S 

1000 
500 



500 1000 1500 2000 

M (GeV) 

Figure 2: The shaded area shows the parameter space of the branon model with a 
thermal relic in the range : tt Br h 2 = 0.122 - 0.098, a contribution to the muon anomalous 
magnetic moment: a M = (29 ± 18) x 10 -10 and favored by the CDMS data. The lower area 
is excluded by single-photon processes at LEP together with monojet signals at Tevatron 
0[5]. Jet and missing energy analysises perfomed with LHC data are the most constraining 
for heavy branons (intermediate area). Prospects for the sensitivity at the LHC for real 
branon production are plotted also for the monojet analysis for a total integral luminosity 
of C = 10 5 and total energy in the center of mass of the collision of yfs = 14 TeV. The 
explicit dependence on the number of branons TV is presented, since all these regions are 
plotted for the extreme values N = 1 and N = 7. 



rect search experiments, if these events are due to branon-nucleus coherent 
interactions. 

Although the present low luminosity of the LHC, present constraints from 
ATLAS and CMS are the most constraining for heavy branons, improving 
previous analyses performed with LEP-II and TEVATRON data. In any 
case, for light branons, these analises are still the most important. The same 
situation is found for standard model processes mediated by virtual branon, 
where current LHC measurements are not competitive yet. 

In addition, there are other signatures which can prove or disprove the 
model. It could be possible to detect branons in DM indirect search exper- 
iments |27]. Two branons may annihilate into ordinary SM matter. Their 
annihilation in places like the galactic halo, the Sun, the Earth, etc. could 
produce cosmic rays to be discriminated through distinctive signatures from 
the background. After annihilation, a cascade process would occur and par- 
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tides such as neutrinos, gamma rays, positrons or antimatter may be de- 
tectable through different experiments such as Atmospheric Cerenkov Tele- 
scopes, Neutrino Telescopes or Satellite Detectors. Work is in progress in 
this direction |28| . 
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